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ABSTRACT 



Aims. We want to improve knowledge of the structure of the inner few AU of the circumstellar disk around the nearby T Tauri 
star TW Hya. Earlier studies have suggested the existence of a large inner hole, possibly caused by interactions with a growing 
protoplanet. 

Methods. We used interferometric observations in the N-band obtained with the MIDI instrument on the Very Large Telescope 
Interferometer, together with 10 /jm spectra recorded by the infrared satellite Spitzer. The fact that we were able to determine N-band 
correlated fluxes and visibilities for this comparatively faint source shows that mid-infrared interferometry can be applied to a large 
number of low-mass young stellar objects. 

Results. The mid-infrared spectra obtained with Spitzer reveal emission lines from HI (6-5), HI (7-6), and [Nell] and show that 
over 90% of the dust we see in this wavelength regime is amorphous. According to the correlated flux measured with MIDI, most 
of the crystalline material is in the inner, unresolved part of the disk, about 1 AU in radius. The visibilities exclude the existence of 
a very large (3 - 4 AU radius) inner hole in the circumstellar disk of TW Hya, which was required in earlier models. We propose 
instead a geometry of the inner disk where an inner hole still exists, but at a much reduced radius, with the transition from zero to full 
disk height between 0.5 and 0.8 AU, and with an optically thin distribution of dust inside. Such a model can comply with SED and 
mid-infrared visibilities, as well as with visibility and extended emission observed in the near-infrared at 2 fim. If a massive planet 
was the reason for this inner hole, as has been speculated, its orbit would have to be closer to the star than 0.3 AU. Alternatively, we 
may be witnessing the end of the accretion phase and an early phase of an inward-out dispersal of the circumstellar disk. 
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1. Introduction 
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The classical T Tauri star TW Hya is the most prominent mem- 
ber of the nearby TW Hya association. With its distance of 
. 51 + 4p c (Mamaiek 20 05j and its estimated age of ~ lOMyr 
'. dHoff et a l. 1998; Web b et al.ll 19991) TW Hya is a Rosetta stone 
' for investigating the transition phase from disks to planets. The 
' obj ect shows clear signs of ongoing accretion, but the rate - 
e.g. lMuzeroUe et all dioOO^ derived ~ 4-lO"'°M0 yr"' -is much 
lower than the typical value of lO '^Mpyr'' found f or IMyr 
old classical T Tauri stars in Taurus jCalvet et al.ll2000h . Images 
taken at various wavelengths from the optical to the millimetre 
regime reveal the presence of a nearly face-on optically thick 
dust disk; see Weinbe rger et al.l ( |2002|) and references therein. 
lApai et all (12004) show the presence of polarised disk emission 
in the K-band as close as (X' 1 (~ 5 AU) to the central star by 
using the polarimetric differential imaging mode of the adaptive 
optics system NACO at the VLT. 

A detailed analysis of the spectral energy distribution (SED) 
led ICalvet et alj ( |2002|) to conclude that the inner disk region 
within 4 AU of the central star is almost devoid of dust, with only 
an optically thin layer of small (sub)micron-sized dust grains re- 
maining. Such a region can be produced by some mechanisms of 



Send offprint requests to: tratzka@aip.de 

* Based on observations with the Very Large Telescope 
Interferometer (VLTI, proposal 075.C-0014). 



dust depletion like photoevaporation, inner planet(s), or signif- 
icant grain coagulation. Indeed, millimetre observations of the 
outer disk p oint to a significant fraction of even centimetre-sized 
dust grains (IWilner et alJl2005h . However, even with HST imag- 
ing and sophisticated AO techniques, it was not possible to char- 
acterise the innermost regions of the disk. iRettig et"aLl (|2004|) 
were able to detect rovibrationally excited CO from the inner 
disk but with no evidence of very hot gas. They conclude that 
TW Hya has cleared its inner disk out to a radial distance of ~ 
0.5 AU and is approaching the end o f its accretio n phase . With 
Keck interferometry measurements, lEisner et"ari (|2006l) could 
spatially resolve the dust continuum emission at a wavelength 
of 2 yum. They concluded that the inner disk consists of optically 
thin, submicron-sized dust extending to within 0.06 AU where it 
may be magnetospherically truncated. 

In this paper we present the first interferometric measure- 
ments of the TW Hya disk in the thermal infrared. The ob- 
servations outlined in Sect. |2] have been performed with the 
'MID-infr ared Interfer ometric instrument' (MIDI) attached to 
the VLTI (iLeinert et al.ll2003allbl) . The data analysis is explained 
in Sect. |3] In Sect. |4] the reduction of spectra obtained with the 
'Spitzer Space Telescope' is described. These spectra are used to 
derive the global dust composition of the mid-infrared emitting 
region (Sect. |5]l that serves as input for our modified Chiang & 
Goldreich model introduced in Sect. |6] By fitting the SED and 
the visibility (Sect. |2ll, the geometrical properties of the circum- 
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Stellar disk can be analysed. The results and their implications 
for our picture of this young stellar object are dicussed in Sect. |8] 
Finally, in Sect. |9] we summarise our findings. 

2. Interferometric observations 

For a spatially and spectrally resolved study of TW Hya and its 
dusty environment in the N-band, we observed this source with 
MIDI. This instrument allows an interfero metric combination of 
the light c oming from a pair of tel escopes jL einert e t al.l2003allbl : 
iMorel et al. 2004; Ratzka 2005; Ratzka et al. 2006). 

The two consecutive measurements analysed and discussed 
in this paper were obtained with the baseline formed by the two 
Unit Telescopes Antu (UTl) and Kueyen (UT2) on 28 May 2005 
within the scope of the MIDI 'Guaranteed Time Observations' 
(GTO), i.e. telescope time dedicated to the MIDI consortium. 
The actual baseline length due to the projection on the sky was 
about 50 m and thus provided a spatial resolution or, more pre- 
cisely, fringe spacing of about 0'.'04 at IQfim. This corresponds 
to 2 AU at the distance of TW Hya. A journal of observations is 
given in Table [T] 

2.1. Observing sequence 

After coarse acquisition by the telescopes, chopped images of 
the source were taken. These images with an angular diameter of 
about 2"allow adjustment of the position of the object to a prede- 
termined pixel on the detector, a 320 x 240 pixel Si:As Impurity 
Band Conduction (IBC) array, in order to maximise the overlap 
of both telescope beams for the subsequent interferometric mea- 
surements. 

Then the beam combiner, which produces two interferomet- 
ric outputs of opposite sign, was put into the optical train. An 
important feature of the measuring process with MIDI is the 
possibility of using a dispersive element. We used a prism with 
a low spectral resolution of AA/A ~ 30, together with a 0'.'52 
wide slit. This spectral resolution not only allowed us to ob- 
tain dispersed fringes, i.e. spectrally resolved interferograms, but 
also increased the coherence length to about + 300 ;um, corre- 
sponding to an interferometric field of view of 1'.'2 for the used 
projected baseline length. Imperfect centring on the white-light 
fringe within a few A thus had no adverse effect on the determi- 
nation of the visibility. 

After finding the location of zero optical path difference 
(OPD) by scanning a few millimetres around the expected point 
of path length equalisation, an interferometric measurement with 
self-fringe-tracking was started. In this mode, the piezo-mounted 
mirrors within MIDI are typically used to scan a range in OPD 
of eight wavelengths {A ~ 10 //m) in steps of 2 fj.m. At each step, 
with the corresponding, fixed OPD, an exposure is taken that 
gives the instantaneous value of the interferometric signal from 
8 to 13yum. After each scan of 80yum, the position of the fringe 
packet in the scan is determined and the VLTI delay lines ad- 
justed in order to recentre the fringe packet for the next scan. The 
scans can either be centred on the white-light fringe or at a fixed 
offset from zero OPD. The latter is advantageous for the coher- 
ent data evaluation (Sect. [3]), which measures the spectral mod- 
ulation of the fringe at all wavelengths simultaneously and thus 
is subject to confusion with background emission at zero OPD. 
The output after a scan is the spectrally resolved temporal fringe 
pattern, which gives the fringe amplitude or correlated flux. The 
scans are repeated in a saw-tooth manner typically 100 - 200 
times to increase the statistical accuracy. No chopping is nec- 
essary for these interferometric measurements, because the two 



resulting output signals can be subtracted from each other, which 
gives an efficient first-order background correction. 

Finally, to derive the (wavelength-dependent) visibilities, 
spectra were recorded on the same pixels of the detector as the 
fringe signal by blocking first the light from one, then from the 
other telescope while chopping the secondary mirror of the UTs 
to remove the thermal background radiation. By definition the 
visibility is obtained as the ratio of correlated flux from the in- 
terferometric measurement and total flux from the photometry 
(Sect. [3]). This gives the raw visibility, still suffering from atmo- 
spheric and instrumental correlation losses. 

The observing sequence for an interferometric measurement 
with MIDI on the VLTI thus produces a diffraction-limited im- 
age of the object within a certain wavelength range defined by 
a filter, a spectrum from 8 to 13 jum, and the correlated fluxes, 
i.e. the visibilities, over the same wavelength range. 

2.2. Calibration 

Since the strength of the fringe amplitude can be reduced by var- 
ious instrumental and atmospheric effects, calibrators of known 
intrinsic visibility have to be observed to determine the instru- 
mental visibility, i.e. the raw visibility that would be observed 
on a perfect point source. The instrumental visibility can be un- 
derstood as the wavelength-dependent transfer function of the 
instrument. Due to the additional correlation losses caused by 
atmospheric turbulence, the calibration measurements should be 
timely and spatially as close as possible to those of the scientific 
target. 

Therefore, we obtained measurements of the calibrator 
HD 95272 immediately before and after the measurements of 
TW Hya by performing the same observation sequence that is 
described in Sect. 12.11 The resulting raw visibilities were then 
corrected for the finite size of the calibrator by assuming a uni- 
form disk with the known stellar diameter. Other calibrators ob- 
served in the same mode on the same baseline during the same 
night (see Table [T]l were used only to check the stability of the 
instrumental visibility (Fig.[T]). This uncertainty also contributes 
to the error bars of the final visibility of TW Hya, but it is smaller 
than the uncertainties in the interferometric measurements of the 
scientific target. 

The interferometric calibrators were taken from the 'MIDI 
Calibrator Catalogue'' of 509 stars that are at least 5 Jy bright 
at 10/im and selected for the absence of circumstellar emission, 
disturbing companions, or strong variability. A subsample^ of 
these stars can also be used for absolute flux calibration. With 
the exception of HD 102839, all calibrators in Table [1] are such 
spectrophotometric standards. However, we restricted ourselves 
to measurements of HD 95272 for determining the spectrum. 



3. MIDI Data reduction and results 

3.1. Visibilities 

For the data reduction of the interferometric measurements, the 
software package MIA+EWS^''^ written in C and IDL was used. 
The con-elated flux F™" (A) and the total flux F™"' (A) as well as 
their ratio, i.e. the wavelength-dependent visibility 



VraJA) 



ptoul (^) ' 



(1) 



http://www.eso.org/~arichich/download/vlticalibs-ws/ 
http://www.eso.org/instruments/midi/tools/spectrophot_std.html 
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Table 1. Journal of observations. 



Date of 


Universal 


Object 




Projected Baseline 


Airmass 


Acquisition 


Interferometric 


Photometric 


Observation 


Time 




[Jy] 


[m] 


[deg] 




Frames 


Frames 


Frames 


28-05-2005 


22:22 - 23:34 


HD 102461 


12.6 


50.4 


24.4 


1.20 


1000x4 ms 


8000x18 ms 


2x4000x18 ms 


29-05-2005 


00:20 - 00:44 


HD 102839 




42.4 


39.9 


1.44 


1000x4 ms 


8000x18 ms 


2x4000x18 ms 


29-05-2005 


01:10-01:41 


HD 95272 


9.5 


55.5 


35.7 


1.20 


4000x4 ms 


8000x18 ms 


2x4000x1 8 ms 


29-05-2005 


01:41 -02:19 


TWHya 




49.8 


41.0 


1.30 


6000x4 ms 


12000x22 ms 


2x6000x22 ms 


29-05-2005 


02:19-02:46 


TWHya 




47.9 


42.4 


1.42 


4000x4 ms 


12000x22 ms 


2x6000x22 ms 


29-05-2005 


02:46 -03:12 


HD 95272 


9.5 


51.6 


37.4 


1.72 


1000x4 ms 


8000x22 ms 


2x4000x22 ms 


29-05-2005 


03:41 - 04:07 


HD 139127 


16.3 


54.3 


27.5 


1.05 


1000x4 ms 


8000x18 ms 


2x4000x1 8 ms 


29-05-2005 


05:17 -05:42 


HD 133774 


11.5 


55.9 


35.3 


1.20 


1000x4 ms 


8000x18 ms 


2x4000x18 ms 



Notes: All observations used the baseline UT1-UT2. The length and the position angle of the projected baseline refer to the time of the 
fringe tracki r ig and the airmass to the time of the photometry. The 10/jm-fluxes of the spectrophotometric calibrators have been taken from 
ICohenetalJ ( 119991) . 




Wovelength [^m] 




10 11 
Wavelength [jtim] 



Fig. 1. Calibrated visibility of TW Hya. The solid curve with er- 
ror bars shows the average of the two independent measurements 
of TW Hya and the associated calibrator HD 95272. The instru- 
mental visibilities derived from the individual observations of 
HD 95272 (dashed) and the other calibrators (dashed-dotted) are 
also shown. A smoothing of 5 pixels in wavelength direction has 
been applied to all measurements. The regions between 9.3;um 
and 9.9 fim and above 12.9 /vm are affected by the atmosphere. 



Fig. 2. Spectrum of TW Hya as determine d by M IDI (solid), 
Spitzer (dashed), and by [Weinberger et al.l (|2002|) with Keck 
(dashed-dotted). The errors of the Spitzer data are negligible and 
those of the Keck spectrum are in general smaller by a factor of 
about two when compared to those of the MIDI spectrum. They 
have not been plotted. A smoothing of 5 pixels in wavelength 
direction has been applied to the MIDI results. 



can be determined with these programs. The first branch of the 
data reduction software called MIA is based on power spectrum 
analysis, i.e. the determination of the spatial modulation of the 
fringe in distinct wavelength bands. The necessary data reduc- 
tion steps are described in a tu toria l on the cited w eb pages, as 
well as in lLeinert etall (|2004 and iRatzkal (l2005h . The second 
branch, called EWS, uses a shift-and-add algorithm in the com- 
plex plane, i.e. averages over suitably modified individual expo- 
sures, to obtain the complex vi sibility A(A) ■ e'^^'^K Thus this is a 
coherent and linear algorithm. iJaffd (12004 and the documenta- 
tion on the web page give further insight into the program. 

In contrast to MIA, which takes advantage of the high corre- 
lated flux level around zero OPD, EWS works best if the OPD 
scan of the interferometric measurements is offset by about 5 A 
from zero OPD for the low spectral resolution of the prism. 
Although we used both methods independently to confirm the 



' http://www.strw.leidenuniv.nl/~koehler/MIA+EWS-Manual 
http://www.mpia-hd.mpg.de/MlDISOFT 



visibility results, we prefer the coherent method that is very 
sensitive in detecting weak fringe signals for our case of a 
faint source. After detailed tests for our case (faint source ob- 
served under good weather conditions), we chose to increase 
the sensitivity further by smoothing in the frame domain with 
a Gaussian width of 15 (parameter gsmooth, set by default to 
4) before determining the group delay. Tests with the tempo- 
ral high-pass smoothing width (parameter smooth) introduced 
to suppress background fluctuations showed no improvements 
when using other values than the default 50. 

The resulting calibrated visibilty, derived by averaging the 
ratios of the raw visibilities of TW Hya and the corresponding 
instrumental visibilities determined by observing HD 95272, is 
shown in Fig. [T] The eiTors are mainly contributed by the un- 
certainties in the interferometric measurements of TW Hya and 
not by those of the transfer function. Also plotted in Fig. [T] are 
the instrumental visibilities determined by the calibrator mea- 
surements obtained on the same baseline in that night. The rms 
scatter of these calibration measurements ranges from ~ 0.02 
between 10 /^m and 11 /im to ~ 0.05 at 8/im and 13;um. The 
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On the other hand, one may think of caHbrating F'^°'^ (A) directly 
from the interferometric measurements instead. When compar- 
ing the two methods, we found in the case of TW Hya that the 
resuhing correlated spectra are very similar, but that the error 
bars are larger when deriving the correlated flux from the al- 
ready calibrated quantities V(A) and F'''"' (/I). These increased 
errors simply reflect that the determination of the correlated flux 
by using Eq. (|2]i requires five independent quantities more than 
the direct method. 

The directly calibrated correlated flux is plotted in Fig. [3] 
The silicate feature looks flatter and broader than in the total 
flux, suggesting we see larger particles here and more crystalline 
dust. A detailed analyis is given in Sect.|5] 



Wovelength [fJ^m] 



Fig. 3. Directly calibrated correlated flux. A smoothing of 5 pix- 
els in wavelength direction has been applied. The regions be- 
tween 9.3 fim and 9.9 //m and above 12.9 fim are affected by the 
atmosphere. 

feature in the wavelength range between 9.3 jjm and 9.9 fim re- 
sults from the ozone band of the earth's atmosphere that leads to 
a reduced count rate at these wavelengths. It is also clearly vis- 
bile in the correlated flux (Fig.[3j. Unfortunately, the extinction 
is very sensitive to the airmass so that a proper normalisation 
of the visibility according to Eq. ([T]i is very problematic in this 
wavelength range. 



3.2. Photometry 

The photometric fluxes 7^'°'^' (X) for the visibility calculations 
are derived by subtracting the background by means of chop- 
ping techniques. However, in particular for faint sources, a sec- 
ond step of background subtraction is required for obtaining a 
reliable spectrophotometry. For this we performed a linear fit 
of the off-source background residuals for each wavelength and 
then subtracted this component from the chopping results. 

The absolute spectra F^^'^ (A) can now be derived from 
^totai (-^-j following the standard procedures applied when re- 
ducing mid-infrared spectra. The different airmasses of the sci- 
entific target and the spectrophotometric calibrator HD 95272 
and the temporal variations in the seeing conditions have to be 
taken carefully into account. For the flux calibratio n we utilised 
the ab solute fluxes for HD 95272 as published by Cohen et"aLl 
(119991) . 

The r esulting spectrum is plo tted in Fig. |2] and the Keck 
data from lWeinberger et al.l (|2002|) are presented for comparison. 
Both spectra agree very well both in shape and in the absolute 
fluxes. Also the Spitzer spectrum of TW Hya (Sect. HJ is in very 
good agreement. By either increasing the Spitzer spectrum or re- 
ducing the MIDI spectrum by 10%, the two spectra fall on top of 
each other. Interestingly, this relative error is what we estimate 
for the absolute calibration error of the Spitzer spectrum. 

3.3. Correlated flux 

The wavelength-dependent correlated flux (A) can now be 
derived by a simple multiplication of the quantities determined 
above: 



= V(A)-F'^^'(A) . 



(2) 



4. Spectrophotometry with the Spitzer satellite 

The presented Spitzer spectra are based on the intermediate 
data products (low-resolution: droopres, high-resolution: rsc) 
processed through the Spitzer data pipeline (ver sion S13.2.0). 
Partially based on the SMART software package jHigdon et al.l 
2004), these intermediate data products were further processed 
using spectral extraction tools developed for the 'Spitzer Legacy 
Science Program s' 'Formation an d Evolution of Planetary 
System s' (FEPS) (iMever et al.ll2006h and 'from Cores to Disks' 
(C2D) dEvans et al.1 1200 3*). For the low-resolution observations 
with the short wavelength module (5 - 14 /vm), the spectra were 
extracted using a 6.0 pixel fixed-width aperture in the spatial di- 
mension. The background was subtracted using associated pairs 
of imaged spectra from the two nodded positions along the slit. 
This also subtracts stray light contamination from the peak-up 
apertures and adjusts pixels with anomalous dark current rela- 
tive to the reference dark frames. Pixels flagged by the 'Spitzer 
Science Center' (SSC) pipeline as 'bad' were replaced with a 
value interpolated from an 8-pixel perimeter surrounding the bad 
pixel. 

For the high-resolution spectra, all 'Basic Calibrated Data' 
(BCDs) from both dither positions were combined within the ex- 
traction algorithm. The extraction was performed by integrating 
over a source profile fit in the cross-dispersion direction along 
the slit. The source profile is a template created from standard- 
star observations (including sky measurements). The width and 
centre of this template were adjusted for TW Hya to encompass 
95% of the observed flux. This fitting process was performed by 
using the highest quality data. Once the source profile was fitted, 
a (uniform) local-sky background level, which is wavelength- 
dependent, was estimated. This method also reduces the effects 
of unidentified bad pixels. 

Both the low- and high-resolution spectra were calibrated 
using a spectral response function derived from IRS ('InfraRed 
Spectrograph') spectra and Cohen or MARCS stellar models for 
a suite of calibrators provided by the SSC. To remove any effect 
of pointing offsets, we matched orders based on the point spread 
function of the IRS instrument, correcting for possible flux 
losses. Offsets can also result in low-level fringing at all wave- 
lengths in the high-resolution spectra. We removed these fringes 
using the irsfringe package developed by Lahuis & Booger^ 
(l2003h . We estimate that the relative flux calibration error across 
orders is 1% and the absolute calibration error 10%. 
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Table 2. Overview of dust species used. For each component its lattice structure (A: amorphous, C: crystalline), chemical compo- 
sition, and shape as well as a reference to the laboratory measurements of the optical constants are specified. 



# Species Structure Chemical Formula Shape Reference 



1 


Amorphous silicate (Olivine stoichiometry) 


A 


MgFeSi04 


Homogeneous 


(1) 


2 


Amorphous silicate (Pyroxene stoichiometry) 


A 


MgFeSizOs 


Homogeneous 


(1) 


3 


Carbons 


A 


C 


Homogeneous 


(2) 


4 


Silica 


A 


Si02 


Inhomogeneous 


(3) 


5 


Forsterite 


C 


Mg2Si04 


Inhomogeneous 


(4) 


6 


Clino Enstatite 


C 


MgSiOa 


Inhomogeneous 


(5) 



References : (1) iDorschner et alj ( Il995h : (2) iPreibisch et aP ( Il993h : (3) iHenning & Mutschkd ( Il997h : (4) IServoin & Pirioul dl973h : (5) 
|JageretalJ ( fl99a) 



1.0 



0.0 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




10 11 
Wavelength [urn] 



Fig. 4. Spitzer IRS observations of TW Hya in the 10 //m wave- 
length region. Shown is the low-resolution IRS spectrum (black) 
between 7 and 13.5/vm with the best-fit dust model described 
in Sect. 15 . II overplotted (grey). As an inset the identically scaled 
high-resolution spectrum is shown between 12 and 13.5 fim. The 
flux-level of both spectra is the same, but the inset has been 
shifted down for clarity. Also indicated in the figure are the po- 
sitions of the main crystalline silicate bands and gas lines. 

Table 3. Results of the spectral decomposition of the low- 
resolution Spitzer spectrum with a continuum temperature 
Tc=1500K and a dust temperature Trf=235K. 

Mass fraction [%] 



Component 


0.1 fim 


1.5 fim 


6.0 fim 


Amorphous Olivine 


66 




2 


Amorphous Pyroxene 


18 


3 




Amorphous Carbon 






4 


Silica 








Forsterite 


3 






Enstatite 


4 


1 





5. The composition of the circumstellar material 

5. 1 . The dust composition 

To characterise the composition of the circumstellar dust we 
analysed the Spitzer low -resolution spectra in the 10 //m wave- 



length region; see Fig. [4] and lUchida et al] ( |2004|) . For this we 
applied a dust model and analysis method for determining the 
physical properties of the dust grains contributing to the ob- 
served emission similar to tho se successfully applied in previ- 
ous studies of the lOyum region dvan Boekel et alj2005l) . Though 
dust grains in protoplanetary disks are most likely highly ir- 
regular aggregates containing many different dust constituents, 
we assume that the aggregates are extremely porous and that 
the individual constituents making up the aggregate may inter- 
act with the radiation field as s eparate entities, as in the cas e 
of interplanetary dust particles jHenning & Stognienko|[T996h . 
Therefore, we assume that the observed emission can be rep- 
resented by the sum of the emission of individual dust species. 
Table |2] summarises the dust species used in our analysis, in- 
cluding all dus t species co mmonly identifie d in pro t oplane - 
tary disks, e.g. ISareent et al. (2006), van Boe kel et alj (l2005h . 
iBouwman et al.1 (1200 ll) " or lBouwman et al.l 020070 . 

Crystalline silicates such as forsterite (Mg2Si04) and en- 
statite (MgSiOs) have many strong and narrow resonances in the 
wavelength range covered by the IRS spectrograph. These rota- 
tional/vibrational bands allow for an accurate determination of 
the chemical composition, grain size, and structure of the crys- 
talline dust species. We find, like in the above-mentioned stud- 
ies, that the pure magnesium end-members of the olivine and 
pyroxene families, together with the assumption of an inhomo- 
geneous grain structure and non-spherical grain shape, give the 
best match to the observed spectral features . Here we use a dis- 
tribution of hollow spheres ( Min et a02005h , which gives an ex- 
cellent match to the observed band positions and shapes. 

For the amorphous iron/magnesium silicates the exact com- 
position and grain shape cannot be constrained as well as for 
the crystalline silicates. Models of the broad amorphous sil- 
icate resonance at ~ lOyum are to some extent degenerate. 
Silicate models with varying magnesium-to-iron ratios and dif- 
ferent grai n shapes/h omogeneity can reproduce the observed 
bands, e.g. lMin et al.l (12007) . To allow for a direct comparison 
with previous studies, we assume here that the amorphous sil- 
icates have an equal magnesium-to-iron atomic ratio and stoi- 
chiometrics consistent with olivine and pyroxene. Further, for 
the amorphous iron-magnesium silicates we assume homoge- 
neous, spherical grains for whic h we calculate the absorptio n co- 
efficients using Mie theory, e.g. lBohren & Huffman! d 1983b . and 
the measured optical constants of these materials (see Table |2] 
for references). 

To take the effect of grain size into account, we calculated 
for each of the different grain species the opacities for three dif- 
ferent grain sizes of 0.1 /im, 1.5//m, and 6.0yum, respectively. 
These sizes sample the range of spectroscopically identifiable 
grain sizes at infrared wavelengths at the SNR of our data. 
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For the analysis of the dust composition we fitted the fol- 
lowing emission model to the Spitzer low-resolution spectra by 
using a linear least-square minimisation; 



5.3. Location of the dust species 



Fy - By(T cont) Co + Bv(rdust) 



3 6 



'J' 



(3) 



where Bv(7^cont) denotes the Planck function at the characteristic 
continuum temperature Tcont, Bv(7^dust) the Planck function at the 
characteristic grain temperature Tdust, Ky is the mass absorption 
coefficient for species j (six in total) and grain size / (three in 
total), and Co and C, j are the weighting factors of the continuum 
and the dust components, respectively. 

Figure H] shows the resulting best-fit model to the Spitzer 
spectrum. As one can see, an excellent agreement between 
model and observations can be obtained using the dust species 
listed in Table |2l without the need for additional components. 
However, the degeneracy between the contribution of the car- 
bons and that of the continuum cannot finally be solved. The re- 
sulting model parameters are listed in Table [3] The mid-infrared 
emission is dominated by small amorphous silicate grains. The 
mass fraction of crystalline grains is about 8%. In the modelling 
of the disk we adopt the dust composition as listed in Table[3j 

5.2. Gas composition 

Apart from the emission features from dust grains, we can also 
observe gas emission lines in the Spitzer spectra (Fig. |4]i. In 
the low-resolution spectra emission lines at 7.47, 12.38 and 
12.81 yum can be detected. We identify these lines with emission 
from HI (6-5), HI (7-6), and [Nell], respectively. The last two 
lines can also be observed in the high-resolution spectrum (see 
inset Fig.|4li. 

iPascucci et al.l ( l2007h observed a number of optically thick 
dust disks in systems having similar accretion rates and X-ray 
luminosities to TW Hya. They detected [Nell] and HI (7-6) to- 
ward one transition disk and [Ne II] alone from three other disks 
that have reduced mid-infrared emission in comparison to clas- 
sical T Tauri star disks. The total line luminosity log L[Neii] 
= 28.25 erg/s measured with Spitzer for TW H ya is similar to 
the va lues observed in the sources studied by iPascucci et alj 
( l2007h . They argue that the [Ne II] emission most likely orig- 
inates from neon atoms at the disk surface layer ionised by 
stellar X-ray or extreme UV photons. The [Nell] line traces a 
gas mass of about 10"^ Mq, which could either be in the disk 
surface layer or the optically thin, depleted inner region (see 
Sect. |6|l of the disk. Unfortunately, the resolution of MIDI is 
too low to detect the line clearly in the correlated spectrum 
and thus to compare its luminosity in the innermost parts of 
the disk with the value derived from the total Spitzer spec- 
trum. However, recent measurements obtained with 'Michelle' 
at Gemini spectrally resolved the [Nell] line and finally re- 
vealed its origin (Pascucci, private communication). On the other 
hand, the observed line luminosities of log Lhi(6-5) =28.87 erg/s 
and log Lhi(7-6) = 28.24 erg/s are too strong to be explained by 
an ionised disk surface or magnetospheric accretion flows. The 
HI lines are most likely originating from an accretion shock 
or th e stellar corona. ROSAT and XMM-Newton observations 
dStelzer & Schmitt 2004; Robrade & Schmitt 2006), which can 
be explained well by X-ray emission from a metal depleted ac- 
cretion shock, seem to support this hypothesis. 



As shown by Ivan Boekel et all (|2004|) . the correlated and uncor- 
rected fluxes provided by MIDI can be used to constrain the 
radial distribution of the emitting silicate grains. When com- 
paring the correlated spectrum (Fig. [3j with the total spectrum 
(Fig.|2]i, it becomes obvious that the overall amorphous silicate 
band is absent in the correlated spectrum. This indicates that the 
region populated by the small amorphous grains can be resolved 
by MIDI and thus is localised at distances larger than a few AU 
from the central star. 

However, a weak remnant of the emission band can be seen 
in the correlated flux. By fitting a first-order polynomial between 
8.5 and 12.5 to the correlated spectrum to estimate the contin- 
uum emission, one is left with the residual emission spectrum 
shown in Fig. |5] When comparing this to the Spitzer spectrum 
after subtracting from the latter the fitted amorphous and contin- 
uum contribution, i.e. taking only the infrared emission from the 
crystalline silicates into account, the match is already very good. 
This shows that the crystalline silicates are located much closer 
to the star than the amorphous species and that the emitting re- 
gion can only be resolved marginally with our interferometric 
observations. The absence of crystalline emission signatures in 
the Spitzer spectra at wavelengths longwards of the N-band is 
another hint of this finding (see Sect. 18.4b . 

Unfortunately, the reduced spectral resolution and the worse 
signal-to-noise ratio of the MIDI spectrum with respect to the 
Spitzer result allow no meaningful spectral decomposition of the 
MIDI residual spectrum. However, a comparison shows that the 
correlated flux contains at least two thirds of the crystalline sil- 
icates we see in the Spitzer spectrum and an enstatite fraction 
that is lower by a factor of two (see Fig.|5]). This finding seems 
to indicate that the enstatite grains have a slightly larger spatial 
extension. But this could also be an artifact of the very simple 
continuum subtraction and the much larger uncertainties in the 
MIDI observations compared to the Spitzer spectrum. 

6. Model 

The stellar parameters of TW Hya are given in Table |4l For 
a more realisti c modelling approach we used a Kurucz model 
(iKuruczl IT992h for a star with 4000 K, a surface gravity of 
logg[cms"'] = 4.5, and solar metallicity at the given distance 
instead of a mere blackbody. Although this template does not 
exactly fit the properties of TW Hya, it is a re asonable approach 
that has also been used by .Eisner et al.l (l2006h . 

6.1. Ttie circumstellar disk 

We describe the disk using a model for passive irradiate d flared 
disks based on the ideas of Chiang & GoldreichI (1199 7!). In this 
picture the dusty part of the disk can be described with two lay- 
ers; a disk interior layer and a disk surface layer. The flared shape 
of the disk ensures that the stellar light can irradiate and heat the 
surface layer to the typical temperature of optically thin dust. 
The absorbed energy is then radiated away by the layer; half of 
this energy is radiated away from the disk and is then observable 
as optically thin warm dust radiation with solid state features in 
emission. The other half is emitted down into the disk interior, 
where it is subsequently re-emitted as optically thick featureless 
emission at wavelengths longer than the emission from the sur- 
face layer. 

Most of th e details of the model u sed here are described 
in the paper by iDullemond et aTl (1200 ih . However, we treat the 
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Fig. 7. Left Panel: Comparison of the observed SED (triangles) and the model fit (black) described in the text. The total SED 
(solid thick) consists of the contribution from the star (dashed-dotted), the inner optically thin disk (solid thin), the midplane of the 
optically thick outer disk (dashed), and its atmosphere (dotted). The additional component in the near-infrared is the contribution 
due to scattering. To demonstrate the influence of the disk mass, the SEDs of disks with a mass of O.I6M0 (black), 0.06 (dark 
grey), and 0.03 Mq (light grey) are given after applying a sedimentation parameter of 0.70 (black), 0.73 (dark grey), and 0.75 (light 
grey) to scale the mid-infrared regime. When in the model fit (thick solid) the sedimentation is no longer taken into account, the 
SED is overestimated in the mid-infrared by a factor of ~ 1.5 (thick dashed), but the visibility reproduces the short wavelengths very 
well (see Fig. [8]). Ri sht Panel: SED lo ngward of lOO um for the four scenarios de scribed above. Measurements: The observational 
data are taken from lMekkadenI (Il998h for the visua l, iRucinski & Krautteij (|1983|) and Sitko et al. (2000) for the near-infrared, this 
paper for the mid-infrared. Weinberg er et alj (l2002h . the IRAS data archive. and lWeintraub et alj (119891) for the sub-millimetre, and 
IWilner et alj ilOOO,) for the millimetre wavelength-regime^. 



inner rim differently. In the iDuIlemond et al.l (1200 lb model the 
rim w as included as a vertical wall. The work of llsella & Nattal 
(I2005h has shown that these rims are not vertical, but strongly 
rounded-off, and recently interferomet ric results could be re- 
produced successfully with this model ( llsella et ani2006') . This 
finding has strong consequences for the inclination dependence 
of the infrared flux and even stronger consequences for the in- 
frared interferometric measurements. We therefore include a sort 
of 'rounded- off' rim in our mod el in a similar way to what is 
described by ICalvet et~aLl (|2002|) : normally the pressure scale 
height of the disk Hp is c alculated self-consistently using the 
IChiang & GoldreichI ( 19971) approach. Close to the inner rim 
however, we artificially reduce Hp smoothly, such that it reaches 
at Rin- The width ^ in which this reduction is done - from no 
reduction at Ri = (1 + ^)Rm to full reduction at Ri^ - sets the 
sharpness of the rim. If ^ = 0.1, the rim is rather sharp, while if 



^ = 1 the transition from no dust to normal flaring is very grad- 
ual. Rather than attempting to calculate ^ self-consistently from 
the physics of possible mechanisms that might smooth-out the 
rim, we simply set ^ and consider it a free parameter. The precise 
way in which we reduce Hp between Ri„ and Ri is also ad-hoc: 
if we call i// the reduction factor, i.e. Hp{R) - t//(R) Hp(R), it is 
defined as 



R - Rin 

R\ — Ru 



(4) 



With the exceptio n of the ISyum measurement, these data are iden 
tical to those used bv lCalvet et 51 12002) and Eisner et al. (2006) for 
their SED. In the Q-band we substitute the result of Jayawardhana et al. 
(1999) by the more recent measurement published in lWeinberger et~ 
I I2002D . 
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Fig. 5. Comparison of the continuum subtracted correlated flux 
(light grey, wavelength range of the ozone band ignored) with 
the residual Spitzer spectrum after subtraction of the best-fit con- 
tinuum and amorphous silicate model (black). Also plotted in 
this figure is a crystalline silicate model of the MIDI continuum- 
subtracted correlated spectrum (dark grey). Like in Fig. |4] the 
positions of the main crystalline silicate bands and the observed 
gas lines have been indicated. 

Table 4. Stellar parameters for TW Hya 

Reference 



spectral type 

Tcff 
mass 

luminosity 

radius 

age 

accretion rate 
extended flux 
distance 
Av 



K7V 
4000 K 
0.6 Mo 
0.19 Lo 
0.91 Ro 
5 - 15Myr 
4xlO-"'Mo/yr 
7% at 2 fim 
51 ±4pc 
0.0 mag 



(1) 
(2) 
(2) 

(3) + parallax 

calculated 

(4) 

(5) 

(6) 

(7) 

(1) 



Refere nces: (1) iRucin ski & Krautter (1983) ; (2) ISchmidt-KaleJ 
JT987); (3) 'Calvet et al. (2002); (4) Uchidaetal.' ("20041); (5) 
Muzerolle et al. (2000); (6) Johns-Krull & Valenti (2001); (7) 
Mamaiekl ( I2005I) 



void 




Fig. 6. Sketch of the geometry of the disk and how the inner dust 
rim is modelled. A vertical cross-cut of the disk is shown. 



modelled this by introducing a reduction parameter that reduces 
the surface height somewhat from the computed value. In our 
model we also allow a (global) parameter that almost directly 
scales the flux from the disk and thus the SED (Fig. |7]), but has 
only very little influence on the brightness distribution and the 
visibility (Fig.|8]l, respectively. 

For the outer dusty part of the disk, we use the mixture of 
dust opacities given in Table [3] derived by modelling the Spitzer 
data (Sect. |5]l. Intrinsically, there is no treatment for scattering 
in the model: only absorption and re-emission is included. The 
total abundance of dust is taken to be 0.01 by mass with respect 
to the gas. 

6.2. Scattering 

Despite the fact that the Chiang-Goldreich model described 
above ignores scattering, this effect is, according to Mie the- 
ory, proportional to /I and hence may play an important role 
at visual and near-infrared wavelengths. For the discussion of 
extended emission in the K-band, we thus include a simple 
model of the additional contribution of scattered light to the 
near-infrared emission (Appendix lAb. 

We find that this contribution is only a small effect (see 
Fig. |7]i and reduces the visibility prediction for TW Hya in the 
K-band by about 0.01 . However, for a fully consistent model, the 
heating losses caused by the scattering would have to be taken 
into account by the Chiang-Goldreich model. 



where q is again a parameter we set freely. Therefore, this ap- 
proach enables us to investigate a variety of modifications of the 
geometry of the rim by just setting the two independent parame- 
ters ^ and q. Since we are mainly intere sted in this ge ometry and 
a direct comparison with the model of ICalvet et al.l (2002) , we 
are confident that this procedure is reasonable, especially when 
taking into account that we have only one independently mea- 
sured visibiliy curve. A pictographic representation of the model 
is given in Fig.|6l 

After predicting the pressure scale height of the disk, the sur- 
face height of the disk can then be computed from the pressure 
scale height when the surface density and dust opacity is known 
and if the assumpti on is made that the d ust and the gas are well 
mixed. As noted bv lChiang et al.l (l200lh dust sedimentation can 
lower the surface height for a given pressure scale height. They 



7. Fitting SED and Visibility 

7. 1 . Previous models 

In the mid-infrared region, where the interferometric measure- 
ments on the VLTI contribute directly, the SED of TW Hya 
(Fig. |7]i is characterised by a strong silicate emission, by quite 
notable excess emission in the range of 3 - 8 um, and by a strong 
increase in Fy towards the 20 - 25 jum range. ICalvet etalJ ( 120021) 
argue that they cannot obtain the strong silicate emission from a 
disk atmosphere, if they do not neglect scattering. They put the 
edge of the optically thick disk at 3 - 4 AU, where it does not 
contribute to the lOfim emission. The silicate emission in their 
model is due to an optically thin distribution of micron-sized 
particles, populating a geometrically rather thin volume in the 
inner disk with constant surface density, which adds sufficiently 
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to the emission longwa rd of the silicate ban d. The 2 fim interfer- 
ometry of TW Hya by lEisner et"^ (l2006l) added the constraint 
that much o f the 7 ± 4% of extended em ission observed at that 
wavelength dJohns-KruU & Valentill200l1) must not come from 
close to the evaporation limit (Rq ~ 0.025 AU), since otherwise 
the visibility could not be as low as their observed value of 0.88 
± 0.05 . They modified the optically thin part of the Calvet et"aD 
(|2002|) model, so that the surface density reflects a radial distri- 
bution according to a power law of 1 .5 and the particles have an 
opacity proportional to /I"' to create both sufficient 2 fim radia- 
tion and the 3-8 yum near-infrared excess. They added a dust-free 
hole inside 0.06 AU, which is required to reproduce the compar- 
atively low K-band visibility. Since they did not try to reproduce 
the SED beyond 8 fim, the lOyum range, where we measured vis- 
ibilities, is not covered by their model. 

Th erefore, we only have to check the model of ICalvet et al.l 
(I2OO2I) for compatibility with the new N-band interferomet- 
ric data. Lacking their original model intensity data, we simu- 
lated them with a Chiang-Goldreich model, keeping position and 
shape for the inner edge of the optically thick disk (3-4 AU), 
constant optical thickness, temperature law, and particle size 
range (0.9 - 2.0 fim) for the optically thin inner disk as given in 
their paper. We assumed the dust composition to be 'astronomi- 
cal silicate' dLaor & Drainel[T993h with a 7% addition of carbon 
for both parts of the disk. The SED predicted from this coarse 
remake looked to be reasonably reproduced and, after adjust- 
ment of the free parameters of our model (see TablejS), fitted the 
observations well, except for wavelengths longer than 100/im. 
We therefore co nsider this s imulation as a valid representation 
of the model of 'Calvet et alj (120021) for the mid-infrared range 
around 10 yum. The visibilities corresponding to this model are 
compared to the observations in Fig. [8] Clearly, the predicted 
values are far too high and show the wrong trend with wave- 
length. To a similar degree the model oy erpredicts the visibil- 
ity observed at 2.14 yum bv lEisner et al.l ('2006). These results 
originate almost exclusively from the physically s imple, easy- 
to-mo del, inner optically thin part of the model of ICalvet et al.l 
( I2OO2I) . which should be well represented by our simulation. In 
spite of its strength in modelling the outer disk, the model will 
apparently need substantial modification to accomodate the new 
interferometric data. Or the other way around, the visibility will 
help to decide whether a model that resembles the SED as well 
as the model introduced bv lCalvet et al.l (|2002|) is more reliable. 

7.2. New mid-infrared constraints 

The new constraints given by our 10 /im interferometry of 
TW Hya will now be used to derive a model of the circum- 
stellar disk geometry for TW Hya and, in this light, to consider 
the uniqueness of the two earlier conclusions about the disk of 
TW Hya - a dust-free zone within 0.06 AU and the limitation of 
the optically thick disk to radii beyond 3 to 4 AU. We will con- 
centrate on the inner 10 AU or so of the circumstellar disk, since 
this is the region to which the data from MIDI refer. We do so by 
modelling the SED and the visibility with the help of the model 
introduced in the previous section. The dust is described by the 
composition derived with an isothermal fit to the silicate feature 
in Sect. |5] Realistic or not, even in our Chiang-Goldreich type 
model with a certain temperature distribution, this composition 
can well reproduce the shape of the silicate emission that is so 
central in our observed wavelength range (see Fig.|9]l. However, 
we are aware of the approximative character of this model. 

Also , we do not want to compete with the detailed argumen- 
tation in lCalvet et al.l (l2002h . which has led to valuable physical 



insights into the outer disk regions. Especially, the explicit dis- 
cussion therein about the disk mass cannot be taken up again, 
because in our dust composition mm- and cm-sized dust grains 
are not considered. These large particles contribute mainly to 
the flux in the sub-millimetre and millimetre wavelength regime 
and thus are not relevant to our mid-infrared analysis. However, 
when including the mm- and cm-sized dust grains, the high 
disk mass of O.I6M0, which is now required to fit the fluxes 
longward of lOOyum, would be reduced significantly. Figure [T] 
demonstrates the ambiguous nature of the mass as a fit parame- 
ter for the mid-infrared when adjusting the sedimentation. The 
differences between the various mass-sedimentation scenarios 
only play a role logward of 100/im. Moreover, the resulting 
mid-infrared visibilities are indistinguishable from one another 
(Fig. [8]). The visibility is thus not able to solve this ambiguity. 
Conclusions based on visibility measurements are basically ge- 
ometric in nature and therefore qualitatively valid and reliable as 
far as questions of disk geometry are concerned. 

In this vein let us investigate the observed visibility curve of 
TW Hya. The in general low values (~ 0.3 from 9 jum to 1 1 yum) 
mean that the typical size of the emitting region - characterised, 
e.g., by the diameter of a structure - is of the order of A/2B 
to A/B, where B is the projected baseline of the interferometric 
measurement, or 1 - 2AU. The increase at short wavelengths 
is a natural consequence of the fact that an increasing fraction 
of the emission comes from unresolved regions, including the 
star. Most informative is the rise of visibility with wavelength 
beyond 1 1 fim. This is the signature of an emitting region whose 
size varies little with wavelength, such as could be produced by 
a 'cliff-like' inner edge of the optically thick outer part of the cir- 
cumstellar disk. Indeed, any disk geometry compatible with the 
visibility measurements seems to require the presence of a 'wall' 
at this inner edge with a diameter of about 1.5 AU and with re- 
duced emission inside. A Chiang-Goldreich-type model of the 
circumstellar disk with an inner radius similar to the sublima- 
tion radius, for example, immediately approaches the observed 
SED, but completely fails to reproduce t he low observed vi sibili- 
ties (see Fig.[3i. The model presented bv lCalvet et al.l (12002 ) also 
fails to reproduce both the shape and the low values of the mea- 
sured visibility, although it includes the required 'cliff-like' rim 
at an even larger distance from the central star than in our model. 
The reason for this finding becomes apparent when analysing the 
contributions of the various disk components to the SED. In the 
model developed by Calvet et al. (2002) the 'wall' is so far away 
from the central heating source that it cannot contribute signif- 
icantly to the N-band flux at all. In this wavelength range the 
visibility traces only the inner disk. 

Although the large fraction of stellar luminosity that is rera- 
diated at longer wavelengths leaves no doubt that the disk ex- 
tending outside this circular wall is optically thick, the proper- 
ties of the region inside this optically thick disk are not well- 
constrained. Nevertheless, if this inne r region is respons ible for 
the observed extended 2 jim emission ( lEisner et al. 2006), it can- 
not be empty. But what we assume there to produce the de- 
sired extended emission is largely arbitrary, as it has been in 
earlier m odels. T he presence of accretion, though at a low level 
(Muzerol le et al.l l2000). was used to justify the presence of op- 
tically thin material in this region showing either a) none (little) 
or b) pronounced silicate emission in the mid-infrared. 

The best combined fit to SED and visibility and our preferred 
solution is a model with an inner disk similar to that described in 
Eisner et al. (2006) and a transition region between the inner and 
the outer disk that is much closer to the star when compared with 
the results of ICalvet et al.l (|2002|) . The parameters of this model 
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Table 5. Parameters for the circumstellar disk of TW Hya 



optically thick outer disk 



inclination 
outer diameter 
surface density 
power law 
sedimentation 
mass 
opacity 



0° (<4°) 
140 AU 
1x10"* g/cm^ 

,.-1.5 

0.7 
0.16 Mo 
see text 



(1) 
(2) 

at 1 AU, fit parameter 
(3) 

fit parameter (see text) 
calculated (see text) 
this paper 





geometry of transition region, 


see Equation (13J( 




0.5 AU 


fit parameter 




0.8 AU 


fit parameter 


£ 


0.6 


calculated 




2.0 


fit parameter 


optically thin inner disk 



^0 


0.06 AU 


fit parameter 


T{Ro) 


1050 K 


fit parameter 


surface density 


9x10"' g/cm- 


at 1 AU, fit parameter 


power law 


,.-1.5 


fit parameter 


mass 


1.2x10^' g 


calculated 


T^(2^ra) 


0.06 


at Rq, calculated 


Ti(10;Um) 


0.013 


at i?o, calculated 


opacity 


10' cmVg 


at 2 yum, (4) 


power law 




(4) 



Re ferences: ( 1) [W einberger et al. (2002); (2) Kri st et alj ( I2OO0I) : 
(3) ID' Alessioetal.1 I I 1 99811 : (4) Eisner et alj ( l200^ 




Wovelength [^m] 



Fig. 8. Comparison of observed visibility and the best simulta- 
neous fit (solid) to both SED and visibility. Without sedimen- 
tation the model overestimates the SED (see Fig. |7]i, but espe- 
cially at short wavelengths the visibility is very well reproduced 
(dashed). Lik ewise, the vis i bihy f ound by reproducing the model 
described by ICalvet et~aLl (l2002i) is shown (dashed-dotted), as 
well as that originating from a pure Chiang-Goldreich model fit- 
ting the SED too (dotted). 



are given in Tables |4]and|5] It is compared to the observations in 
Figs.|2]and[8l 



8. Discussion 

8.1. The inner disk 

We have no good reasoning for the character of the inner disk. 
We only know that its 3 - 8 fim emission must be substantially 
reduced with respect to a standard Chiang-Goldreich flaring disk 
model. 

Between the two possibilities considered by us, we ob- 
tain the best representation of the data by assuming an opti- 
cally Jhindistrihutiotj of particles showing the opacity given in 
Eisn er et'al] ^2006), i.e. proportional to /I"'. In analogy to so- 
lar system interplanetary dust dynamics, we estimate that the 
Poynting-Robertson lifetimes even for 1 mm particles at the tran- 
sition region would be only 10^ yr dWvatt & Whipplelll950h . For 
this model we therefore have to assume that accretion continues 
through the inner edge of the outer disk and that (smaller) parti- 
cles showing features in the mid-infrared are removed on a much 
shorter timescale by a combination of typical sources also act- 
ing in interplanetary space: radiation and stellar wind pressure, 
magnetic acceleration, collisions, gas drag, and the Poynting- 
Robertson effect. All of these are more important for small par- 
ticles. 

Our second scenario is an optically thin distribution of sili- 
cate dust showing pronounced silicate emission. This does not 
produce a very good fit to the SED and visibility data because it 
tends to overpredict the lOyum emission. Moreover, the timescale 
for the Poynting-Robertson drag is proportional to the radius of 
the particle, which only makes this scenario feasible if there is 
continuous replenishment of these grains on short timescales. 

However, both scenarios would face the same problem of re- 
producing the low observed K-band visibilities if the inner disk 
continues inward to the dust sublimation limit at ~ 0.025 AU 
(~ 1500 K), because a too large fraction of the 2/im radiation 
would then be unresolved by the 62 m basel ine, i.e. u = 49.6 m 
and V = 36.7 m, of the Keck observat ions (lEisner e t al. 200^. 
Therefore, one is led to postulate, like lEisner et al.l (i2006) did, 
a truncation of the inner disk at Rq ~ 0.06 AU or ~ 12 stel- 
lar radii, probably caused by the action of the stellar mag- 
netic field, but somewhat out side the corotation radius of 6.3 R* 
(iJohns-Krull & ValentilllOOTl) . Knowledge of the inner boundary 
of the inner disk may profit in the same way from spectrally dis- 
persed near-infrared interferometry as the knowledge of the in- 
ner edge of the outer disk benefits from the interferometric MIDI 
observations in the mid-infrared. This complementary observa- 
tion could already be perf ormed with the sen sitivity provided 
by AMBER on the VLTI (IPetrov et al.l l2003l) or the success- 
fully tested grism mode of the Keck Interferometer jEisner et al.l 
2007). 

Qualitati vely we thus reproduce the model results of 
lEisner et al.l (l2006l) . Quantitatively, we have to assume a particle 
density that is higher by a factor of ~ 1.5 to have the same in- 
frared excess. However, our prediction for the squared visibility 
at a waveleng th of 2. 14 u m an d for the projected baseline length 
reported in E isner et al. (12006) is 0.92 and thus consistent with 
the measurements at the Keck interferometer. 

The amount of diffuse 2 jim emission of r — 0.07 + 0.04 

dJoims-Krull & Valenti 2001) can be kept in both scenarios dis- 
cussed above, but the first one - comparatively large particles 
- allows the best overall representation. We therefore propose 
to consider this model for the inner disk as a serious possi- 
bility that deserves further investigation. We refer to the work 
by Tanaka et al. (2005) and note that Cissla (2007) discussed a 
mechanism that can lead to strong particle growth in the inner 
part of a circumstellar disk. 
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8.2. Transition region 

The visibility measurements at 10 fim put the edge of the outer 
disk at 0.5 - 0.8 AU. This is much less than the radius of 3 - 4 AU 
for the inner disk derived by Calvet et al. (2002). However, it re- 
mains the question what interaction caus es this abrup t change 
between the outer and inner disk. iCalvet et al.l (|2002|) present 
arguments that a massive planet might be the underlying rea- 
son. They also mention GG Tau and DQ Tau as examples of 
where a companion has cleared the inner region of a larger disk 
- documented by millimetre- and near-infrared-images of the 
GG Tau system - and where accretion from th e outside still oc- 
curs - at least in the case of DQ Tau. Although lWeinberger et al.l 
(1200 2) searched for such a companion to TW Hya and excluded 
any point source brighter than a planet of lOMj for separations 
> 50 AU from TW Hya, this result is not directly relevant for 
our model, where t he orbit o f a planet or companion would have 
to be at ~ 0.25 AU. lHerbig' & Bell ( 1988) noted that three radial 
velocity determinations showed a spread of +11 km/s, making 
it worthwile to check for a close companion. A companion of 
0.1 Mg at 0.20 - 0.25 AU from TW Hya would have an orbital 
velocity of 43 - 48 km/s leading to a radial velocity of the pri- 
mary of 0.5 - 0.6 km/s when assuming an orbital inclination of 
4° (see Table Is). This is much less than the scatter in these early 
observations. We are not aware of radial velocity monitoring of 
TW Hya, although such a search could be fruitful. The low near- 
infrared brightness of a 0. 1 Mq companion will make it difficult, 
but not impossible, to detect it eve n by an interferom etric instru- 
ment like AMBER on the VLTI dPetrov et al.ll2003l) . However, 
we emphasise that locally enhanced particle growth or photoe- 
vaporation of the dust can also lead to an optically thin inner 
disk. 

According to our analysis, the abrupt change from the 
outer to the inner disk, which was the main feature found by 
ICalvet et all ( 12002 ). occurs approximately by a factor of six 
closer to the star than previously concluded. The 10 //m visibil- 
ities put hard constraints on the location of the inner edge of 
the outer disk. It happens that this new, smaller distance essen- 
tially coincides with the region of CO emission at about 0.5 AU. 
This region was identified by 'Rettig et al] (|2004|) on the basis 
of the observed rotational temperature of CO transitions. The 
edge of the outer disk, if located at the small radius derived by 
us, would n a turally qualify as the origin of the CO emission. 
ICalvet et al.l (l2002h proposed the large radius of the optically 
thin inner disk due to their conclusion that the inner disk and 
not the upper layers of the optically thick outer disk produce 
the strong s i licate emission feature. However, this solution of 
ICalvet et al.l (|2002|) is not unique but just one of the possibilities 
for describing the SED unconstrained by mid-infrared visibili- 
ties and the mm-images. We showed in Sect. |7.2| that, fully con- 
sistent with the mid-infrared visibilities, the silicate emission can 
be produced in the atmosphere of the outer optically thick disk 
when taking moderate sedimentation into account. 

8.3. Disk evolution 

With its estimated age of ~ lOMyr, TW Hya is among the old- 
est classical T Tauri stars. We might expect that its circumstel- 
lar disk shows signs of evolution, like continuing crystallisation, 
growth of particles, beginning dispersal of the disk, and onset of 
planet formation. 

The evidence is mixed. Given the l arge size of the c ircum- 
stellar disk visible in scattered light (IWeinberger et al.l |2002|) 
and the strongly disk-dominated emission for wavelengths above 



10 fim, which led lCalvet et alJ (l2002l) to a model with a disk mass 
of 0.06 Mq, the disk around TW Hya is far from dispersal, even 
after lOMyr. Calvet et al. (2002) provide interesting suggestions 
to explain this unusual survival time. On the other hand, the 
existence of mm- to cm-sized particles , inferred from the slow 
fall of mm-emissio n with wavelength dWeinberger et al.l 120021 : 
IWilner et al ] 120051) . shows strong particle growth into a size- 
regime where dust settling to the midplane of the disk should 
become effective. 

If most of the 10 yum emission comes from the upper disk 
layers, as we have argued, a much more primordial size distribu- 
tion is present there. The mid-infrared measurements in general 
would not be sensitive to mm-sized particles, but the lOfim sil- 
icate emission feature would become broader and flatter in the 
prese nce of micron-sized par ticles. This occurs in many T Tauri 
stars ("Przvg odda et alj 12003!) . but not in TW Hya. Instead the 
emission profile here requires a vastly dominating fraction of 
submicron-sized silicate particles. Since general particle growth 
by coagulation would have not only led to an increase in the 
fraction of large grains in the midplane, but also in the disk at- 
mosphere, our analysis indicates processes, e.g. sedimentation, 
that have removed these particles. 

A second mode of dust evolution is crystallisation. Again, 
very little crystalline material contributes to the observed sili- 
cate emission from TW Hya. In the interferometrically observed 
flux, representing emission from the innermost region with a di- 
ameter of a few AU, and probably coming from the inner edge of 
the disk, the situation is different. This is not totally unexpected 
since a similar behavio ur has been found in Herbig Ae/Be stars 
dvan Boekel et alJ l2004). but it is not obvious that the same has 
to hold for the much less luminous T Tauri stars. 

From these observations, it appears that both radial mixing 
that would distribute crystals and vertical mixing that would 
homogenise size distributions are not important in the disk of 
TW Hya. For a deeper discussion of transport mechanisms it will 
be helpful to have the corresponding data available for a number 
of T Tauri stars. 



8.4. Location of tiie crystalline grains 

A direct indicator of an inhomogeneous distribution of the crys- 
talline grains in the optically thick disk of TW Hya are the sili- 
cate emission features at wavelengths longer than 20 yum. These 
features are present in the model SED (Fig.]?!, but they are ab- 
sent in the Spitzer spectrum. When removing from the entire 
disk the crystalline species, one finds the measured smooth SED 
at longer wavelengths, but the N-band that causes the dust com- 
position of the model is no longer well represented (see Fig.|9]l. 
Interestingly, when keeping the found dust mixture with its low 
percentage of crystalline grains only in the transition region, 
both the N-band emission and the featureless SED at longer 
wavelengths is reasonably reproduced. 

Although there is evidence - both from the correlated flux 
and the SED - that the processed dust is close to the edge of 
the optically thick disk and may even be distributed throughout 
the optically thin inner disk, detailed investigations of the dust 
composition as a function of distance from the star are necessary 
for deeper insight into the location of the crystalline grains. We 
leave it to other studies to conclude on this. 
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Fig. 9. Mid-infrared part of the SED for three different distribu- 
tions of the crystalline grains: from the dust mixture found by 
fitting the Spitzer spectrum the crystalline grains have been re- 
moved nowhere (solid), outside the transition region (dashed), 
and in the entire disk (dotted). The underlying triangles repre- 
sent the measured fluxes (for references see caption of Fig. |7]l. 

8.5. Millimetre maps 

The images obtained at 7 mm bv IWilner etaP d200(]h are not 
discussed here. They ar e well reproduced b y the outermost part 
of the disk modelled by lCalvet et al.l (12002 *). which in principle 
could be joined to the modified inner parts of the disk, whose 
properties have been determined in this work with the help of 
mid-infrared measurements. 

9. Summary and conclusion 

To summarise and conclude: 

- We were able to interferometrically measure the compara- 
tively faint T Tauri star TW Hya. Similar studies will there- 
fore be possible for a large number of low-mass young stars. 

- We confirm that the circumstellar disk of TW Hya can be 
described best by a rather standard outer part with an abrupt 
transition to an inner part with reduced emission. However, 
based on lOyum visibility measurements, we put this transi- 
tion region much closer to the star than it has been done in 
earlier models, namely between 0.5 and 0.8 AU. This transi- 
tion region would qualify as the source of much of the CO 
emission observed by .Rettig et al.. (2004) to originate from a 
region around 0.5 AU. 

- The inner disk, as already discussed bv lEisner et alj (l2006l) . 
does not appear to extend inwards fully to the dust sublima- 
tion and corotation limit. Spectrally dispersed, near-infrared 
interferometry may be needed to definitely clarify this point. 

- Our results do not contribute anything new to the ques- 
tion of what causes the abrupt change in properties between 
the outer and the inner part of the disk. Possible answers 
are changing opacities induced by locally enhanced particle 
growth or a massive planet. However, in our proposed geom- 
etry, any potential companion or massive planet would have 
to revolve closer to the star and faster, and should thus be 
easier to detect by radial velocity monitoring. 

- Comparison of spectrophotometry from Spitzer with spec- 
trally resolved correlated flux measurements from MIDI at 



the VLTl shows that most crystalline material appears to 
be concentrated within a few AU from the star. The disk of 
TW Hya is not mixed well at the present epoch. 
- The parametrised models used in this work and in earlier 
papers are not completely satisfactory for application to an 
even modestly complex disk like the one of TW Hya. A fully 
fledged radiative transfer/hydrodynamics code might give a 
physically more convincing representation of the now copi- 
ously available data relating to the geometry of this circum- 
stellar disk. 
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Appendix A: Scattering 

We calculate the contribution of light scattered at the optically 
thick part of the disk by 

Isca, U) = F, (A) T (A) A (A) O {A, m , (A. 1) 

where is the irradiation by the star at the location where the 
scattering occurs, t ~ 1 the optical thickness of the scatter- 
ing surface, A the albedo, an d O the phase function taken from 
iHenvev & Greensteinl(ll94ll) as 

O {A, J?) = h _^ ^ (^^2 _ 2g (^) cos M'^'^ (A.2) 

4 ;r L J 

with the average cosine of the angle of scattering g and the scat- 
tering angle § ~ 90°. The values for A and g are those for a 
mixture of uncoate d submicron-sized silicate and graphite par- 
ticles published bv lMathis et al.l (Il983l) . The foreground extinc- 
tion has been assumed to be zero (see Tabled at all wavelengths. 
This is suflicient for the near-infrared, but might be erroneous for 
shorter wavelengths. We further restricted the diameter of the 
emitting region to five times the Airy-disk of a 10 m telescope at 
the corresponding wavelength. 
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